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Abstract 

Recent tests at Tyndall Air Force Research Laboratory show that an elastomer- 
coated to the inside wall of a building can offer significant protection for 
occupants by keeping wall fragments together and shielding out blast debris. In 
this paper, we develop an equivalent single degree-of-freedom model that can be 
used to predict dynamic response of a polymer-retrofitted concrete brick wall 
subjected to a stand-off explosion. As an example, we consider the blast 
response of a 10 fi-square, 8 in-thick concrete brick wall coated with a 0.083 in¬ 
thick layer of polyurea. Explosive field tests with a uniformly-distributed 
pressure pulse of 40 psi peak value and 20 ms pulse duration indicate that the 
retrofitted wall reaches a maximum deflection of about 7 in. Our analytical 
model suggested that the blast response of this wall could be simulated assuming 
the concrete modulus was reduced to about 5% of its original value. The 
proposed single degree-of-freedom model, which is based on coupling of the 
bending/membrane resistance of the wall, compared very well with the 
ABAQL1S results when the maximum deflections of the wall were between 1-2 
times the wall thickness. 

Keywords: blast response, elastomer-coated concrete wall, equivalent single 
degree-of-freedom model. 

1 Introduction 

A major cause of injury and death during bomb attacks on buildings is flying 
debris due to fragmentation of walls, windows, furniture and equipment. The 
United States Air Force Research Laboratory at Tyndall Air Force Base recently 
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conducted full-scale tests showing that concrete masonry unit (CMU) walls and 
trailers retrofitted with polyurea or polyurethane coatings will protect occupants 
when these structures are subjected to a standoff explosion [1], By bonding 
tightly to the wall and being able to undergo large membrane deflections without 
tearing, these elastomers mitigate damaging effects of the blast. The cracked 
wall is kept together by the elastomer and shields against flying debris from the 
outside. The ability of the elastomer to stretch extensively without breaking is 
extremely important because it allows the blast energy to be transferred into 
viscoelastic strain energy of the wall. This blast energy would have otherwise 
been used in fragmentation damage - turning pieces of wall, window, furniture 
and equipment into deadly projectiles. 

In this paper, we present an equivalent single degree-of-freedom model that 
can be used to predict dynamic response of a retrofitted concrete brick wall 
subjected to a stand-off explosion. The use of equivalent or lumped parameter 
single degree-of-freedom systems to obtain the blast response of structures is 
covered extensively in Biggs [2]. This practical approach also forms the basis of 
blast design codes in the United States [3]. Fragmentation damage of the 
concrete bricks under blast is a very complex phenomenon. Not only is the 
initial deformation and fracture properties of concrete sensitive to strain rate and 
pressure [4], but comminution and post-rubblization flow affect damage 
evolution in a concrete wall. An average concrete modulus equal to some 
fraction of the actual modulus of concrete is thus assumed. The reduced concrete 
modulus is a parameter that can be determined either empirically or with more 
refined FEA and rate-dependent continuum damage model for concrete [5]. 

2 Problem formulation 


Consider a 2ax2a elastomer-coated CMU wall that is clamped at the base and 
pin-supported at the top. This wall is subjected to a stand-off explosion, which 
we will describe as a uniformly distributed pressure pulse on the wall: 


P(t) = P 0 


h-T 

V 'ey 


0 ) 


where p 0 is the peak force and x is the decay time of the wall. The retrofitted 

wall is modeled by the equivalent single degree-of-freedom (SDOF) system 
shown in Fig. 1, where A is the deflection at the center of the wall, 

F eq = 4a p(t) is the equivalent force of the blast, M eq is the lumped mass of the 

CMU and elastomer, P eq is a function describing the bending and membrane 

resistance of the wall, and the dots denote derivative with respect to time. The 
equation of motion for this SDOF system is 

M eq A + P eq (A)=F eq (t) (2) 
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Pinned 

a x edge 



Fully clamped Equivalent System 

edge 


Figure 1: Single degree of freedom model elastomer coated CMU wall 

subjected to blast. 


Lumped or equivalent mass and spring forces are calculated by assuming the 
following displacement and velocity fields: 


w(x) = A 1 


-V 

UJ 


(3) 


w(x) = A 


1 - 






(4) 


where the dot denotes time derivative and the origin of the x-y coordinate system 
is in the center of the wall. 


2.1 Equivalent mass 


Equating the kinetic energy of an equivalent mass M eq to the total distributed 
kinetic energy of wall gives 


M ec, =4a(p CMU h, +p c hJ p V - .~ dx =^-a 2 (p CMU h l +p e h e ) 

J A“ to 


32 

15 


(5) 


where p CMU is the average mass density of the CMU block (including any 
hollow sections), p t , is the density of the elastomer, and h^and h e are the 
thickness of the CMU block and elastomer, respectively. 
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2.2 Bending-membrane resistance function 


The bending-membrane resistance function or equivalent spring force P cq 

shown in Fig. 1, is found using the principle of minimum potential energy to 
derive an expression for the static resistance of the wall subjected to a uniformly 
distributed pressure p 0 . The elastic strain energy of the wall is 


a 

U = 2a J(N x e xo +M x K x )dx 


( 6 ) 


where N x and M x are the effective membrane force and bending moment of the 


composite wall, respectively, e xo =- 


dw 

dx 


is the midline axial strain for 


moderately large deformations and k x = 


d 2 w 

dx 2 


is the curvature. The effective 


membrane force and bending moment of the composite wall depend on the 
modulus of the concrete and elastomer and the geometry of the CMU block. 
Consider a standard CMU block with outer dimensions h,x h( xb,and two 

hollow sections of dimensions h 2 x h 2 x b 2 , as shown in Fig. 2. The distance to 
the neutral axis of the composite wall, h NA , is located by setting the resultant 
axial forces equal to zero: 


J 


h N'A 

jVxb +<**m )dzdy = 0 
—(hi h na ) 


(7) 


where the bending a xb and membrane a xm stresses are 

“(h i+h e -h NA )<z<h NA -h e 

h na ~ h e <z<h NA 

"(hi +h e -h NA )< z < h NA -h e 

h NA -h e < z < h NA 

and E c and v c are the Young’s modulus and Poisson’s ratio of the cement used 
to make the CMU block, respectively; E e and v e are the Young’s modulus and 

Poisson’s ratio of the elastomer, respectively; and the z-axis is defined with 
respect to the neutral axis, as shown in Fig. 2. 
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Figure 2: Cross-sectional geometry of CMU block and elastomer. 


Substituting cr xb and cr xm into Eq. (7) and integrating gives the following 

expression for the position of the neutral axis in combined bending and 
membrane tension: 



An effective bending moment for the composite system is found by 

b,/2 h NA 

M x = — j* |a xb zdzdy = Dk x (9) 

1 —b, / 2 —(h, +h c -h NA ) 
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The effective membrane force is also found by 

. b, n h NA 


N x=t- J Ja xm dzdy = Ce X( 


( 10 ) 


1 -b, / 2 — (h, ► li c Iina ) 

(b,h, — 2b 2 h) E c h £ 
where C = ^-4- ~—^-r — r~T l' 

b, m) M) 

Substituting Eqs. (9) and (10) into (6) and using Eq. (3) to integrate the resulting 
function over x, gives the following strain energy: 


U=- 


■ (d,a 2 + d 2 a’ + D 3 A 


3 . r, a 4 1 , 8C a4 


\ 8C aH 

) + 7T A 

5a 


where 


Dl - L 2 \[( h NA0 — ^ e ) _ (hl +h e _ ^ NAO ) ] 

i l i ~ V c / 


2b 2 E c a 


h NAO ~ 


r h,~h^ 


A 3 


" h e 

J J 


+ 


3b,(l-v c . 

77" ~yt[3h e h NA0 -3h“h NA0 +h e 
3 v — v c ; 


Y h, +h 2 ^ 

A 2 “ 


A 3 


+ h e h NA() 


01 ) 


D, =■ 


E c a 

(Ml 


2b 2 E c a 

3b,M) 


(hNAO h e ) (h,+h e h NA0 ) 


hNAO “ 

V V 


h, -h 2 


A A 

2 

f h 1 + h 2 'l 

. A 

1 

01 

“ h e 

— 

1 l 

+ h e ~h NA() 


/ ) 


l 2 J 

9 



+ 


E e a L , 

j [2h nao h e 


3(l-V, 2 ) 
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D, = 


5(1 -v c 2 ) b, 5(1 -vj)' 

The potential of the work done by the external loads is 

V = ~fjp° wdxdy = ~~7 ~ A 


( 12 ) 
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The total potential energy is 


n = u + v=- 


•(d,a 2 + d ; 


A 



(13) 


Minimizing the total potential energy with respect to A gives the following 
pressure equilibrium equation: 




10D!A + 15D 2 A 2 +4(aC + 5D 3 )A 3 ] (14) 


Thus the equivalent force equilibrium equation is 


P eq =Po( 2 a ) 2 =-^-[lOD,A + 15D 2 A 2 + 4(aC + 5D 3 )a 3 ] (15) 

5a 


3 An example 


As an example, consider the blast response of a 10 ft-square wall constructed 
from 8x8x16 standard CMU blocks [6] and coated on the distal side of the 
explosion with a 0.083 in-thick layer of polyurea. Explosive field tests with a 
blast pressure of peak value p 0 = 40 psi and pulse duration of t = 20 ms indicate 
that this wall will reach a maximum deflection of about 7 in. 


Geometric parameters of the composite wall are as follows: 


a = 60 in 
h, = 7.625 in 
h 2 =5.125 in 
b, = 15.625 in 
b 2 =6.3125 in 
h e = 0.083 in 


half-width and half-height of wall 

outer block height 

height of block hollow section 

outer block base 

base of block hollow section 

elastomer thickness 


Mechanical properties of the concrete are as follows: 

E c = 3,000 ksi Young’s modulus of concrete 

v c =0.15 Poisson’s ratio of concrete 

Pc mu = 7.2x10 5 lbs 2 /in 4 mass density of homogenized CMU block 
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Mechanical properties of the polyurea are taken from experiments [7]: 


fEejS 

0<£<8y 


|(E e i -E e2 )e y + E e2 e 

e> e y 


[l7.5 + 27.86e 

0 < s < 7 

[3 + 4.058 

where E pl = < 

61 [212.5 

, E e - 

>={ 

8>7 

■ [10 


0 < e < 1.7 
e> 1.7 


f0.16-0.042s 

jo.oi 


ins 


-l 


0 < s < 3.6 

, and E e . 

e >3.6 


and E e2 are in ksi and e is 


v e = 0.5 Poisson’s ratio of elastomer 

p e = 9.076xl0 -5 lbs 2 /in 4 mass density of polyurea 

Substituting the equivalent mass, equivalent spring force, and equivalent 
loads into Eq. (2) gives a nonlinear, nonhomogeneous second-order differential 
equation, where C, D],D 2 and D 3 are rate-dependent functions. This 

differential equation is solved with initial conditions A = 0 and A = 0 and 
reduced concrete modulus in the range of 7.25-145 ksi, using MATLAB. These 
reduced concrete moduli were chosen so that the maximum deflection of the wall 
would be roughly between 1-2 times the thickness of the wall since this is the 
range in which the bending-membrane tension would occur in the wall. At 
E c =145 ksi, the maximum deflection of the wall is about 7 in, which is the 
maximum deflection that was obtained from the blast test. Since the original 
modulus of concrete is E c = 3,000 ksi, cracking of the concrete in the CMU 
block must have caused the modulus to decrease by about 4.8% of its original 
value during the blast. A parametric study on the influence of the elastomer 
modulus indicates that wall deflections are fairly insensitive to elastic properties 
of the polyurea. Thus the only important role of the elastomer is to keep 
concrete fragments together so that the blast energy can be absorbed in bending 
and membrane stretching of a “weakened” concrete wall. Without the elastomer, 
CMU blocks would have separated into flying debris and the wall would have no 
resistance to the blast. 

Since the response of the wall is not strongly dependent on mechanical 
properties of the elastomer, we assume the concrete and the elastomer as two 
rate-independent, linear elastic materials and use ABAQUS Explicit to find the 
central deflection of the elastomer-coated wall in our example problem. The 
modulus of the elastomer is chosen as E e =32.5 ksi and all other material 

parameters except for the reduced concrete modulus are kept the same as in the 
example. Variations of the maximum deflection with the reduced modulus of 
concrete from ABAQUS Explicit and the SDOF model are compared in Fig. 3. 
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The analytical model over-predicts FEA results by about 58% when the 
maximum deflection from ABAQUS is 4.5 in, but the over-prediction decreases 
to about 3.2% when the deflection is 16.5 in. In general, FEA and analytical 
predictions approach each other as the deflection increases from about 8 in, 
which is the thickness of the wall. This is to be expected since the SDOF was 
based on the bending-membrane resistance of the wall and this theory breaks 
down if deflections are too small (less than the thickness of wall) or too large 
(many times the wall thickness). 



Figure 3: Comparison of predicted maximum deflection from equivalent 

single degree-of-freedom (SDOF) model and ABAQUS Explicit 
(E e =32.5ksi). 

4 Conclusions 

An equivalent single degree-of-freedom model of an elastomer-coated wall 
subjected to blast was presented in this paper. Equivalent mass, force resistance 
and loads were derived by assuming parabolic shape functions for the transverse 
deflection and velocity of the wall. As an example, the problem of a blast on a 
10 fit-square wall constructed from 8x8x16 standard CMU blocks and coated on 
the distal side of the explosion with a 0.083 in-thick layer of polyurea, was 
considered. Explosive field tests with a blast pressure of peak value 40 psi and 
pulse duration of 20 ms indicated that this wall reaches a maximum deflection of 
about 7 in. It was found that the concrete modulus must have been reduced to 
145 ksi or about 5% of its original value during the blast. 
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The analytical model over-predicted FEA results by about 58% when the 
maximum wall defection was about half of the thickness, but it was only about 
3.2% higher than the ABAQUS solution when the deflection was about twice the 
wall thickness. Since the proposed single degree-of-freedom model was based 
on coupling of the bending and membrane resistance of the wall, it should only 
be used when the maximum deflection of the wall is expected to be greater than 
the wall thickness. 
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